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In this paper, we give a proof of the existence of stationary dark soliton solutions or heteroclinic 
orbits of nonhnear equations of Schrodinger type with periodic inhomogeneous nonhnearity. The 
result is illustrated with examples of dark solitons for cubic and photorefractive nonlinearities. 



I. INTRODUCTION 

T— I \ 

, Nonlinear Schrodinger (NLS) equations have a great number of applications in different fields of modern physics, 
from laser wavepackets propagating in nonlinear materials to matter waves in Bose-Einstein condensates (BEC), 
, gravitational models for quantum mechanics, plasma physics or vifave propagation in biological and geological systems 
a : ^' -2-4, 28]. 

^ ^ \ From the physical point of view two of the most relevant forms of the nonlinear equations of Schrodinger type are 
the stationary NLS equation with cubic nonlinearity (namely cubic NLS), 

00 ■ 

^ ; -u" + Xu + au^ ^0 (1) 

' and the stationary NLS equation with cubic-quintic nonlinearity (namely cubic-quintic NLS), 

n, '. 

. ■ -u" + Xu + au^ + hv!^ = 0, (2) 

1—1 where a, 6 e M. In the case when a < 0, for the first model (respectively, 6 < for the second model), and with 

1 ^ I boundary conditions 

lim u{x) = 0, lim u' {x) = 0, (3) 

|2:|— >oo |a;|— >oo 



' the solutions of Eq. ([1]) are called solitary waves or bright solitons. When a > 0, (respectively, 5 > 0) nontrivial 
. dynamics requires that the solution does not vanish at infinity. In this case the so-called dark solitons or heteroclinic 
' orbits appear. 

In the last years there has been an increasing interest in the one-dimensional nonlinear Schrodinger equation with 
a spatially modulated nonlinearity: 

- u" + \u + a{x)u^ = (4) 

^ with a; G M and a{x) describing the spatial modulation of the nonlinearity. This equation arises in different physical 

contexts such as nonlinear optics and Bose-Einstein condensation, being in the later field which has gained more 
interest in the last year. This fact has been triggered by the possibility of u sing the Feschbach resonance management 
5^ , techniques to modify spatially the collisional interactions between atoms [2^ . Different aspects of the dynamics 

of solitons in these frameworks have been considered such as the emission of solitons [2l| and the soliton propagation 
when the space modulation of the nonlinearity is either a random fll , periodic j29j , linear [25| or localized function 

In recent years, the study of existence of bright solitary waves has gained a lot of interest from the mathematical 
point of view (see e.g. 0, [3, H^, H^). The existence of dark solitons has been previously considered in non-periodic 
potentials (see Jjy and references therein). To the best of our knowledge, the only proof of existence in NLS models 
with periodic potentials (optical lattices) has been presented for the cubic-quintic nonlinearity in [13 . 

Many NLS models, including Eq. can be considered as specific cases of a more general equation with a periodic 
nonlinearity coefficient 

-Uxx + )^u + g{x)f[u)u^G, (5) 
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with 5 : R — ^ M being T-periodic, satisfying the fohowing properties: 

< 5min ^ g{x) ^ (6a) 
gi^) = 9{-x). (6b) 

and / : R — > R satisfying 

/eCi(R), /(-w) = /(w) > 0,u e R, and /'(u) > 0, u e R+. (7) 
Usuahy a dark sohton is defined as a solution of Eq. ([S]) verifying the following asymptotic boundary conditions 

u{x) 



u±{x) 



1, X ^ ±oo, (8) 



where the functions u±(x) are sign definite, T-periodic, real solutions of Eq. ([5]). By symmetry of Eq. ([S]), U-{x) = 
—u+{x), which allows to consider odd solutions for u{x) and to reduce complexity of the problem. This kind of 
solutions are also denoted as asymptotically oscillatory kinks or black solitons and possess a topological charge. 
Notice that these solutions are different to bubbles, for which lima;_j.+oo u{x) = lima;_j._oo u{x) [3-1^. 

The aim of this paper is, therefore, to extend earlier studies to the model given by (O, which, to the best of our 
knowledge, has not previously been considered. From the mathematical point of view, the strategy of proof consists 
of several techniques from the classical theory of ODE's (upper and lower solutions) and planar homcomorphisms 
(topological degree and free homeomorphisms) combined in a novel way. 

The rest of the paper is organized as follows. In Section 2 the preliminaries of the paper are introduced and a proof 
of the existence of periodic solutions for Eq. ([5]) is given. Section 3 contains the main result about the existence of 
heteroclinic orbits or dark solitons. In section 4 the theorems of Section 3 are applied to two typical examples: the 
cubic and photorefractive nonlinearities. We conclude in Section 5 with a summary of the main findings of the paper. 

II. EXISTENCE OF PERIODIC SOLUTIONS 

We start by analyzing the range of values of A for which the existence of nontrivial solutions of Eq. ([5]) can be 
found. 

Theorem 1 //A 0, the only hounded solution of Eq. 0^ is the trivial one, u — Q. 

Proof: Let u be a nontrivial solution of Eq. ([5]). Multiplying ([5]) by u and integrating over the whole space, one gets 

{{u:^f + Xu^ + g{x)f{u)u^) dx^O (9) 

Using (|6ap . (|6b[) and, by hypothesis, the fact that A ^ 0, the only bounded solution is u = 0. ■ 

Therefore, throughout this paper, we will take A < 0. As gmin ^ gi^) ^ ffmax, let us consider two auxiliary 
autonomous equations: 

iW+Au(i)+.g„,in/(Ji))7.« =0 (10) 

422+^„(2)+^_^^(j2)y2)^0 (11) 

These equations possess two nontrivial equilibria 

e^ = .r'(-—), (12) 



g 



mill 



i''' = r'{-^). (13) 

We need to classify these points. Linearizing Eqs. ([TU)) and ([TT]) around of equilibria points (|12l) and one obtains 
the following eigenvalues 



- ±J.W(/-(-^))/-(-^), (14) 

.ymin / / \ .ymi 



M^^) = ±W5nW ( f-' ( f-'i^] (15) 
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As / is a continuous and increasing function, there exists f~^. Moreover is also an increasing function. Thus, the 
points ((H]), are saddle points. Moreover, due to the evenness of the function /, there exist two extra equilibrium 
points more, = — C*-^' and = — We would also like to note that S,'^^^ > since that / is increasing on 
the positive semi-axis. 

Before continuing, let us recall a useful definition from the theory of ODE (see, for example Let the following 

second order differential equation be 

Uxx^h{x,u), (16) 
where h is a continuous function with respect to both arguments and T-periodic in x. 
Definition 1 (i) We say that u : [a, +oo) — > K is a lower solution of J J 6)] if 

Uxx > h{x,u) (17) 

for all X > a. 

(a) Similarly, u : [a, +oo) — > M is an upper solution of \lb]) provided that 

Uxx < (18) 

for all X > a. 

We shall now prove the existence of a T-periodic and an unstable solution between both points (^^^^^ and 

Proposition 1 The points and which were previously calculated, are respectively constant upper and lower 
solutions of Eq (0). Moreover there exists an unstable periodic solution situated between them. 

Proof By using Eq. (|10[). we obtain: 

+A^'^^ > \e'^+g^^nf{e^)e'^=Q (19) 

and, similarly, for the Eq. (fTTj) : 

-eE^ + AC^^) +g(:r)/(e(2))e^ < Ae(^) +g™../(^(2))e(^) = (20) 

Thus, by the latter definition, ^^^^ and ^^^^ are a couple of well-ordered upper and lower solutions, respectively. 
Therefore, there exists a T-periodic solution between them [Tl|. Such solution is unstable because the associated 
Brouwer index to the Poincare map is —1 (see, for example [l8|). ■ 

We have therefore a positive and T-periodic solution of Eq. ([5]), u+{x), satisfying ^^^^ ^ w+(-^) =^ ^^^^ ■ Due to the 
symmetry of the equation we also have a negative solution u^{x) = —Uj^{x). 



III. EXISTENCE OF A DARK SOLITON. 



In this section, we prove the existence of a heteroclinic orbit connecting the periodic solutions u_ and u+. This 
heteroclinic orbit can also be called a dark soliton. 

The following theorem, proved in [23] with the use of ideas from [l(| , is the key one for this work. 

Theorem 2 Let w,v : [a, +oa) R be a bounded functions verifying the following conditions 

1. w{x) < v{x), yx > a 

2. Wxx{x) > h{x,w) and Vxx{x) < h{x,v), Va; > a. 
Then there exists a solution u{x) of il6]) such that 

w{x) < u{x) < v{x) (21) 

// moreover, there exists x such that 

3. min > 0, 
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then there exists an T -periodic solution p{x) such that 

lim {\u{x)- p{x)\ + \u.,{x)- p,,{x)\)=Q (22) 



Moreover, p{x) is the unique T-periodic solution in the interval [infa;^a;„ sup^^^.^^ . 

Let us apply this theorem to our model. We have that 

h{x,u) = \u{x) + g{x)f{u{x))u. (23) 

Since g{x) is symmetric, one can consider x ^ and extend the obtained solution u{x) as an odd function for x < 0. 
The solutions of Eqs PH)) and dTTl) . u^-^'' and u'-^\ which are heteroclinic orbits joining — C'^-* with ^'^-'^^ and — with 
respectively, satisfy the conditions ^ and ([2]) of Theorem[2j with v{x) = m*^^-'(x) and w{x) = u'^'^\x). We thus 
have a bounded solution u{x) of Eq. ([5]) such that 

u(2)(a;) < u(a;) < u(i)(x) (24) 

Hence, in order to prove that there exists a solution u{x) of ([5]) converging to u+(x) and U-{x), which were obtained 
in Proposition [TJ as a; — >■ ±oo, one has to verify the condition ([3]) of Theorem [51 
As a can be taken arbitrarily large, condition ([3]) is equivalent to 

min [X + g{x)f{u)+g{x)f{u)u]>Q (25) 

Here, we must distinguish four different cases: 
(i) f'{u) is an increasing function, for u > 0. Then (1251) is now equivalent to 

A + 9minf{e^) + .9mi„/'(e^'^)e^'^ > (26) 

Taking into account Eqs. (IT^ and p^ . a straightforward analysis of the last inequality gives the following 
estimate for A 

A - ( r\-^)] r\-^) > (27) 



As / is an increasing function, condition (j27p can be satisfied for certain values of gmin, 5max, A and, therefore, 
it is a sufficient condition for the existence of dark solitons. 

(ii) f'{u) is a decreasing function, for u > 0. Then, constraint ([25]) is reduced to 

A + 5mi„/(^(')) + gnunf'i^^'^)^^'^ > (28) 

and subsequently to the following inequality 

A - A^ + .9„,„/' f /-i(-^)) r\-^) > (29) 

Again, condition (j29p can be satisfied, depending on the values of ffmin, ffmax and A. 

(iii) f'{u) has, at least, a minimal point in ug G C*'^^]- For this case, the constraint (^5)) is reduced to 

A + .9mi„/(e^'^) + .9mi„/'("o)e^'^ > (30) 

and hence 

A-A^+5mi„/'(uo)r'( ^)>0 (31) 

5'max 

Thus ([5T|) is a sufficient condition for the existence of dark soHtons. 
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(iv) f'{u) has a maximal point in ui € 

[^(2)^^(1)]. Thus, (HH) is reduced to 

A + gmin/(e^'^) + <?mi„/'(C^'^)e^'^ > (32) 

if we supposse that /'(^''^-') < /'(C'^"*)- Otherwise, we hold 

A + 5mi„/(^('>) + gnunne'^)^'^ > (33) 

Thus, Eq. p2p gives the following estimate for A 

A - A^ + ,g,„i„/' f/-i(-^)) f-\-^) > (34) 
and, alternatively, Eq. ([55]) gives the following estimate for A 

A - A^ + g„,„/' (r'i- — )) f-'h^) > 

^max \ ^min / ^max 

which are Eqs. (H?]) and (^5]) . 



(35) 



IV. APPLICATION TO THE EXISTENCE OF DARK SOLITONS 

In this section we check the consequences of our main result in the problem of the existence of dark solitons to the 
stationary nonlinear Schrodinger equation with periodic nonlinearity, given by Eq. ([5]) . 

We distinguish hereafter two cases of great physics relevance: the NLS with either cubic or photorefractive nonlin- 
earity. 

A. Dark solitons in the cubic NLS equation 

The cubic nonlinearity, f{u) — u^, is by far the most common in several fields of physics, as nonlinear optics 
[isj (where it is also dubbed as Kerr nonlinearity) and Bose-Einstein condensation [l^. Dark solitons arise in self- 
defocusing optical Kerr media [l^ or repulsive Bose-Einstein condensates [l^ . 

Thus, if f{u) = V?, Eq. (O can be cast as 

— Uxx + Am -I- g{x)u'^ — (36) 
Therefore, we are in the case (i) of the previous section. Thus, condition Q is equivalent to 

min [2A + Qg(x)v^] > (37) 

This last inequality is equivalent to 2A-|-6gmin('?^^^)^ > 0. Using Eq. (IT^ and the fact that A < 0, we find a connection 
between ^min and gmax: 

5min > — (38) 

Consequently, if (j38p is verified, dark solitons will exist in the cubic nonlinear Schodinger equation with periodic 
nonlinearity. 

We consider below an example of dark (black) soliton of the Eq. ([55)) . in order to check the concepts introduced in 
Section 3. This example corresponds to the periodic nonlinearity g{x) given by 

gi^) = — ^ (39) 

^ ' (1 + Q;cos(a;a;))3 ^ ' 

with u] ~ 2-y/[Af, 50 and a < 1 being positive constants. In order to satisfy the relation psp . a must fulfill the 
constrain a < (3^^'^ — l)/(3^/'' + 1). For this case, the solutions of Eqs. (|TU)) and (|TT|) . for f{u) = v? can be found 
analitically, being 



= J— tanhl , /^tanhlJ^x) (40) 

V .9min \ V 2 / y gmax \ V 2 / 
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FIG. 1: [Color Online]. The dark soliton u{x) (solid blue line), the hyperbolic periodic solutions u±{x) (dashed red and green 
line), and the heteroclinic orbits u^-^^ and it'^' (dashed black line) connecting ±^'^' and ±^'^' respectively, for the parameters 
A = —0.5, go — 1, a — 0.1 of the cubic NLS equation. 



Following 0, [1] , the solution of Eq. (15^ with the boundary conditions (|^ is 



, , IjJ / 1 — Q;2 _ 

u(x) = —\ V 1 + acoswxtanh 

2 V 50 



uj 1 — 



2 V 2 



X{x) 



(41) 



with X{x) being a solution of 



tan(-Vl-a2x(.)j^^^tan- 



(42) 



Then, the boundary condition u+ is 



(jj j [1 — a'^){l + acosLox) 



90 



(43) 



and u_ = — u+. 

This solution is depicted in the Fig. [TJ 



B. Dark solitons in the photorefractive NLS equation 



There exists many optical materials, like lithium niobate, that exhibits photorefractivity. This nonlinear effect, 
contrary to the Kerr media, implies saturation, i.e. the function / is increasing and fulfills [23|: 

/(0) = 0, lim /(u) = /oo<oo (44) 

The study of photorefractive materials is of high interest, as possess unique features. For instance, they are able to 
exhibit either self-focusing or self-defocusing in the same crystal; besides, solitons can be generated with low optical 
power p^.[l7l|. 
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FIG. 2: [Color Online] . The dark soliton u{x) calculated for the parameters A = —0.3, Qmin ~ 2, pmax = 8 of the photorefractive 
NLS equation. 



Photorefractive nonlinearity is given by 



and, consequently, Eq. ^ transforms into 



/(") = Y 



Uxx + Au + g{x)— — 5- = 
1 + 



Then, for / given by Eq. (j45p . we have one maximum and one minimum for f'{u), given respectively by 



Thus, we are considering the case (iii) of the previous section. Then, Eq. (j3ip becomes 



A - A 1- — — 



mm 

5'max o V ''^ ymax 



A + gn 



> 



If we expand expression (1481) . we obtain the following inequality 



A^ + 5maxA 



27 (ffminffn 



which can be written as 



64 (5min - ffmax)^ 

(A-Ai)(A-A2) >0 



> 



(45) 



(46) 



(47) 



(48) 



(49) 



(50) 
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where 



^ ffmax ^ ffmax 27 5min 



16 (5min - 5max)^ 



^ ffmax _ ffmax /^^ 27 ffmin 



provided that 



> 1 



16 (^min ^max)'^ 



27 



(51) 
(52) 

(53) 



Thus, from ([SO)), we obtain that at least for 



Ae (-ffmax, A2] U [Ai,0) 



(54) 



Eq. (|46)) has dark sohton solutions. 

Analytical solution of Eq. (|46l) cannot be generally found. A numerical solution for g{x) = 2 + 6cos^(x) and 
A = —0.3 is depicted in Fig. [21 For this case, gmin = 2, (/max = 8. We want to note that for these values, the inequality 
(1551) is satisfied, being Ai = —0.39 and A2 = —7.60. Thus, the sufficient condition ([Ml) is satisfied as A belongs to the 
interval given by (1541) . 

In order to find the numerical solution, it has been implemented a fixed point Newton-Raphson algorithm with a 
finite-differences discretization of the Laplacian, supplemented by anti-periodic boundary conditions. 



V. CONCLUSIONS 



In this paper, we have considered the existence of dark solitons or heteroclinic orbits of the nonlinear Schrodinger 
equation with periodic nonlinearity. The proof method relies on some results of the qualitative theory of ordinary 
differential equations that require some concepts such as upper and lower solutions, topological degree or free home- 
omorphisms. As an example, we have constructed an analytical black soliton-solution for cubic nonlinearity and a 
numerical black soliton for photorefractive nonlinearity. 
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